home *** CD-ROM | disk | FTP | other *** search
/ Enter 2004 January / enter-2004-01.iso / files / maxima-5.9.0.exe / {app} / share / maxima / 5.9.0 / src / numerical / slatec / zunhj.lisp < prev    next >
Encoding:
Text File  |  2003-02-09  |  60.8 KB  |  1,072 lines

  1. ;;; Compiled by f2cl version 2.0 beta 2002-05-06
  2. ;;; 
  3. ;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
  4. ;;;           (:coerce-assigns :as-needed) (:array-type ':simple-array)
  5. ;;;           (:array-slicing nil) (:declare-common nil)
  6. ;;;           (:float-format double-float))
  7.  
  8. (in-package "SLATEC")
  9.  
  10.  
  11. (let ((ar (make-array 14 :element-type 'double-float))
  12.       (br (make-array 14 :element-type 'double-float))
  13.       (c (make-array 105 :element-type 'double-float))
  14.       (alfa (make-array 180 :element-type 'double-float))
  15.       (beta (make-array 210 :element-type 'double-float))
  16.       (gama (make-array 30 :element-type 'double-float))
  17.       (ex1 0.33333333333333337)
  18.       (ex2 0.6666666666666667)
  19.       (hpi 1.5707963267948966)
  20.       (gpi 3.141592653589793)
  21.       (thpi 4.71238898038469)
  22.       (zeror 0.0)
  23.       (zeroi 0.0)
  24.       (coner 1.0)
  25.       (conei 0.0))
  26.   (declare (type double-float conei coner zeroi zeror thpi gpi hpi ex2 ex1)
  27.            (type (simple-array double-float (30)) gama)
  28.            (type (simple-array double-float (210)) beta)
  29.            (type (simple-array double-float (180)) alfa)
  30.            (type (simple-array double-float (105)) c)
  31.            (type (simple-array double-float (14)) br ar))
  32.   (f2cl-lib:fset (f2cl-lib:fref ar (1) ((1 14))) 1.0)
  33.   (f2cl-lib:fset (f2cl-lib:fref ar (2) ((1 14))) 0.10416666666666669)
  34.   (f2cl-lib:fset (f2cl-lib:fref ar (3) ((1 14))) 0.08355034722222221)
  35.   (f2cl-lib:fset (f2cl-lib:fref ar (4) ((1 14))) 0.12822657455632716)
  36.   (f2cl-lib:fset (f2cl-lib:fref ar (5) ((1 14))) 0.29184902646414046)
  37.   (f2cl-lib:fset (f2cl-lib:fref ar (6) ((1 14))) 0.8816272674437576)
  38.   (f2cl-lib:fset (f2cl-lib:fref ar (7) ((1 14))) 3.3214082818627677)
  39.   (f2cl-lib:fset (f2cl-lib:fref ar (8) ((1 14))) 14.995762986862555)
  40.   (f2cl-lib:fset (f2cl-lib:fref ar (9) ((1 14))) 78.92301301158652)
  41.   (f2cl-lib:fset (f2cl-lib:fref ar (10) ((1 14))) 474.4515388682643)
  42.   (f2cl-lib:fset (f2cl-lib:fref ar (11) ((1 14))) 3207.490090890662)
  43.   (f2cl-lib:fset (f2cl-lib:fref ar (12) ((1 14))) 24086.549640874004)
  44.   (f2cl-lib:fset (f2cl-lib:fref ar (13) ((1 14))) 198923.1191695098)
  45.   (f2cl-lib:fset (f2cl-lib:fref ar (14) ((1 14))) 1791902.0077753437)
  46.   (f2cl-lib:fset (f2cl-lib:fref br (1) ((1 14))) 1.0)
  47.   (f2cl-lib:fset (f2cl-lib:fref br (2) ((1 14))) -0.14583333333333334)
  48.   (f2cl-lib:fset (f2cl-lib:fref br (3) ((1 14))) -0.09874131944444445)
  49.   (f2cl-lib:fset (f2cl-lib:fref br (4) ((1 14))) -0.14331205391589505)
  50.   (f2cl-lib:fset (f2cl-lib:fref br (5) ((1 14))) -0.31722720267841353)
  51.   (f2cl-lib:fset (f2cl-lib:fref br (6) ((1 14))) -0.9424291479571203)
  52.   (f2cl-lib:fset (f2cl-lib:fref br (7) ((1 14))) -3.5112030408263544)
  53.   (f2cl-lib:fset (f2cl-lib:fref br (8) ((1 14))) -15.727263620368046)
  54.   (f2cl-lib:fset (f2cl-lib:fref br (9) ((1 14))) -82.28143909718595)
  55.   (f2cl-lib:fset (f2cl-lib:fref br (10) ((1 14))) -492.35537052367056)
  56.   (f2cl-lib:fset (f2cl-lib:fref br (11) ((1 14))) -3316.218568547972)
  57.   (f2cl-lib:fset (f2cl-lib:fref br (12) ((1 14))) -24827.67424520859)
  58.   (f2cl-lib:fset (f2cl-lib:fref br (13) ((1 14))) -204526.58731512981)
  59.   (f2cl-lib:fset (f2cl-lib:fref br (14) ((1 14))) -1838444.91706821)
  60.   (f2cl-lib:fset (f2cl-lib:fref c (1) ((1 105))) 1.0)
  61.   (f2cl-lib:fset (f2cl-lib:fref c (2) ((1 105))) -0.20833333333333337)
  62.   (f2cl-lib:fset (f2cl-lib:fref c (3) ((1 105))) 0.125)
  63.   (f2cl-lib:fset (f2cl-lib:fref c (4) ((1 105))) 0.3342013888888889)
  64.   (f2cl-lib:fset (f2cl-lib:fref c (5) ((1 105))) -0.40104166666666674)
  65.   (f2cl-lib:fset (f2cl-lib:fref c (6) ((1 105))) 0.0703125)
  66.   (f2cl-lib:fset (f2cl-lib:fref c (7) ((1 105))) -1.0258125964506173)
  67.   (f2cl-lib:fset (f2cl-lib:fref c (8) ((1 105))) 1.8464626736111112)
  68.   (f2cl-lib:fset (f2cl-lib:fref c (9) ((1 105))) -0.8912109375)
  69.   (f2cl-lib:fset (f2cl-lib:fref c (10) ((1 105))) 0.0732421875)
  70.   (f2cl-lib:fset (f2cl-lib:fref c (11) ((1 105))) 4.669584423426247)
  71.   (f2cl-lib:fset (f2cl-lib:fref c (12) ((1 105))) -11.207002616222994)
  72.   (f2cl-lib:fset (f2cl-lib:fref c (13) ((1 105))) 8.78912353515625)
  73.   (f2cl-lib:fset (f2cl-lib:fref c (14) ((1 105))) -2.3640869140625)
  74.   (f2cl-lib:fset (f2cl-lib:fref c (15) ((1 105))) 0.112152099609375)
  75.   (f2cl-lib:fset (f2cl-lib:fref c (16) ((1 105))) -28.212072558200244)
  76.   (f2cl-lib:fset (f2cl-lib:fref c (17) ((1 105))) 84.63621767460074)
  77.   (f2cl-lib:fset (f2cl-lib:fref c (18) ((1 105))) -91.81824154324002)
  78.   (f2cl-lib:fset (f2cl-lib:fref c (19) ((1 105))) 42.53499874538846)
  79.   (f2cl-lib:fset (f2cl-lib:fref c (20) ((1 105))) -7.368794359479632)
  80.   (f2cl-lib:fset (f2cl-lib:fref c (21) ((1 105))) 0.22710800170898437)
  81.   (f2cl-lib:fset (f2cl-lib:fref c (22) ((1 105))) 212.57013003921713)
  82.   (f2cl-lib:fset (f2cl-lib:fref c (23) ((1 105))) -765.2524681411817)
  83.   (f2cl-lib:fset (f2cl-lib:fref c (24) ((1 105))) 1059.9904525279999)
  84.   (f2cl-lib:fset (f2cl-lib:fref c (25) ((1 105))) -699.5796273761325)
  85.   (f2cl-lib:fset (f2cl-lib:fref c (26) ((1 105))) 218.1905117442116)
  86.   (f2cl-lib:fset (f2cl-lib:fref c (27) ((1 105))) -26.491430486951554)
  87.   (f2cl-lib:fset (f2cl-lib:fref c (28) ((1 105))) 0.5725014209747314)
  88.   (f2cl-lib:fset (f2cl-lib:fref c (29) ((1 105))) -1919.457662318407)
  89.   (f2cl-lib:fset (f2cl-lib:fref c (30) ((1 105))) 8061.722181737309)
  90.   (f2cl-lib:fset (f2cl-lib:fref c (31) ((1 105))) -13586.550006434138)
  91.   (f2cl-lib:fset (f2cl-lib:fref c (32) ((1 105))) 11655.393336864532)
  92.   (f2cl-lib:fset (f2cl-lib:fref c (33) ((1 105))) -5305.646978613403)
  93.   (f2cl-lib:fset (f2cl-lib:fref c (34) ((1 105))) 1200.9029132163525)
  94.   (f2cl-lib:fset (f2cl-lib:fref c (35) ((1 105))) -108.09091978839464)
  95.   (f2cl-lib:fset (f2cl-lib:fref c (36) ((1 105))) 1.7277275025844574)
  96.   (f2cl-lib:fset (f2cl-lib:fref c (37) ((1 105))) 20204.291330966145)
  97.   (f2cl-lib:fset (f2cl-lib:fref c (38) ((1 105))) -96980.5983886375)
  98.   (f2cl-lib:fset (f2cl-lib:fref c (39) ((1 105))) 192547.00123253153)
  99.   (f2cl-lib:fset (f2cl-lib:fref c (40) ((1 105))) -203400.17728041552)
  100.   (f2cl-lib:fset (f2cl-lib:fref c (41) ((1 105))) 122200.46498301746)
  101.   (f2cl-lib:fset (f2cl-lib:fref c (42) ((1 105))) -41192.65496889756)
  102.   (f2cl-lib:fset (f2cl-lib:fref c (43) ((1 105))) 7109.514302489364)
  103.   (f2cl-lib:fset (f2cl-lib:fref c (44) ((1 105))) -493.915304773088)
  104.   (f2cl-lib:fset (f2cl-lib:fref c (45) ((1 105))) 6.074042001273483)
  105.   (f2cl-lib:fset (f2cl-lib:fref c (46) ((1 105))) -242919.1879005513)
  106.   (f2cl-lib:fset (f2cl-lib:fref c (47) ((1 105))) 1311763.6146629772)
  107.   (f2cl-lib:fset (f2cl-lib:fref c (48) ((1 105))) -2998015.9185381066)
  108.   (f2cl-lib:fset (f2cl-lib:fref c (49) ((1 105))) 3763271.297656404)
  109.   (f2cl-lib:fset (f2cl-lib:fref c (50) ((1 105))) -2813563.226586534)
  110.   (f2cl-lib:fset (f2cl-lib:fref c (51) ((1 105))) 1268365.2733216248)
  111.   (f2cl-lib:fset (f2cl-lib:fref c (52) ((1 105))) -331645.1724845636)
  112.   (f2cl-lib:fset (f2cl-lib:fref c (53) ((1 105))) 45218.76898136273)
  113.   (f2cl-lib:fset (f2cl-lib:fref c (54) ((1 105))) -2499.8304818112097)
  114.   (f2cl-lib:fset (f2cl-lib:fref c (55) ((1 105))) 24.380529699556064)
  115.   (f2cl-lib:fset (f2cl-lib:fref c (56) ((1 105))) 3284469.853072038)
  116.   (f2cl-lib:fset (f2cl-lib:fref c (57) ((1 105))) -1.9706819118432228e+7)
  117.   (f2cl-lib:fset (f2cl-lib:fref c (58) ((1 105))) 5.095260249266464e+7)
  118.   (f2cl-lib:fset (f2cl-lib:fref c (59) ((1 105))) -7.410514821153266e+7)
  119.   (f2cl-lib:fset (f2cl-lib:fref c (60) ((1 105))) 6.634451227472903e+7)
  120.   (f2cl-lib:fset (f2cl-lib:fref c (61) ((1 105))) -3.756717666076335e+7)
  121.   (f2cl-lib:fset (f2cl-lib:fref c (62) ((1 105))) 1.3288767166421817e+7)
  122.   (f2cl-lib:fset (f2cl-lib:fref c (63) ((1 105))) -2785618.1280864547)
  123.   (f2cl-lib:fset (f2cl-lib:fref c (64) ((1 105))) 308186.4046126624)
  124.   (f2cl-lib:fset (f2cl-lib:fref c (65) ((1 105))) -13886.089753717042)
  125.   (f2cl-lib:fset (f2cl-lib:fref c (66) ((1 105))) 110.01714026924674)
  126.   (f2cl-lib:fset (f2cl-lib:fref c (67) ((1 105))) -4.932925366450996e+7)
  127.   (f2cl-lib:fset (f2cl-lib:fref c (68) ((1 105))) 3.2557307418576575e+8)
  128.   (f2cl-lib:fset (f2cl-lib:fref c (69) ((1 105))) -9.394623596815784e+8)
  129.   (f2cl-lib:fset (f2cl-lib:fref c (70) ((1 105))) 1.55359689957058e+9)
  130.   (f2cl-lib:fset (f2cl-lib:fref c (71) ((1 105))) -1.6210805521083374e+9)
  131.   (f2cl-lib:fset (f2cl-lib:fref c (72) ((1 105))) 1.1068428168230145e+9)
  132.   (f2cl-lib:fset (f2cl-lib:fref c (73) ((1 105))) -4.958897842750303e+8)
  133.   (f2cl-lib:fset (f2cl-lib:fref c (74) ((1 105))) 1.4206290779753308e+8)
  134.   (f2cl-lib:fset (f2cl-lib:fref c (75) ((1 105))) -2.4474062725738724e+7)
  135.   (f2cl-lib:fset (f2cl-lib:fref c (76) ((1 105))) 2243768.1779224495)
  136.   (f2cl-lib:fset (f2cl-lib:fref c (77) ((1 105))) -84005.4336030241)
  137.   (f2cl-lib:fset (f2cl-lib:fref c (78) ((1 105))) 551.3358961220206)
  138.   (f2cl-lib:fset (f2cl-lib:fref c (79) ((1 105))) 8.147890961183122e+8)
  139.   (f2cl-lib:fset (f2cl-lib:fref c (80) ((1 105))) -5.866481492051847e+9)
  140.   (f2cl-lib:fset (f2cl-lib:fref c (81) ((1 105))) 1.8688207509295826e+10)
  141.   (f2cl-lib:fset (f2cl-lib:fref c (82) ((1 105))) -3.4632043388158773e+10)
  142.   (f2cl-lib:fset (f2cl-lib:fref c (83) ((1 105))) 4.1280185579753975e+10)
  143.   (f2cl-lib:fset (f2cl-lib:fref c (84) ((1 105))) -3.302659974980072e+10)
  144.   (f2cl-lib:fset (f2cl-lib:fref c (85) ((1 105))) 1.79542137311556e+10)
  145.   (f2cl-lib:fset (f2cl-lib:fref c (86) ((1 105))) -6.563293792619285e+9)
  146.   (f2cl-lib:fset (f2cl-lib:fref c (87) ((1 105))) 1.5592798648792575e+9)
  147.   (f2cl-lib:fset (f2cl-lib:fref c (88) ((1 105))) -2.251056618894153e+8)
  148.   (f2cl-lib:fset (f2cl-lib:fref c (89) ((1 105))) 1.7395107553978165e+7)
  149.   (f2cl-lib:fset (f2cl-lib:fref c (90) ((1 105))) -549842.3275722887)
  150.   (f2cl-lib:fset (f2cl-lib:fref c (91) ((1 105))) 3038.0905109223845)
  151.   (f2cl-lib:fset (f2cl-lib:fref c (92) ((1 105))) -1.4679261247695616e+10)
  152.   (f2cl-lib:fset (f2cl-lib:fref c (93) ((1 105))) 1.1449823773202579e+11)
  153.   (f2cl-lib:fset (f2cl-lib:fref c (94) ((1 105))) -3.990961752244665e+11)
  154.   (f2cl-lib:fset (f2cl-lib:fref c (95) ((1 105))) 8.192186695485775e+11)
  155.   (f2cl-lib:fset (f2cl-lib:fref c (96) ((1 105))) -1.0983751560812234e+12)
  156.   (f2cl-lib:fset (f2cl-lib:fref c (97) ((1 105))) 1.008158106865382e+12)
  157.   (f2cl-lib:fset (f2cl-lib:fref c (98) ((1 105))) -6.453648692453764e+11)
  158.   (f2cl-lib:fset (f2cl-lib:fref c (99) ((1 105))) 2.879006499061506e+11)
  159.   (f2cl-lib:fset (f2cl-lib:fref c (100) ((1 105))) -8.786707217802326e+10)
  160.   (f2cl-lib:fset (f2cl-lib:fref c (101) ((1 105))) 1.763473060683497e+10)
  161.   (f2cl-lib:fset (f2cl-lib:fref c (102) ((1 105))) -2.167164983223795e+9)
  162.   (f2cl-lib:fset (f2cl-lib:fref c (103) ((1 105))) 1.4315787671888897e+8)
  163.   (f2cl-lib:fset (f2cl-lib:fref c (104) ((1 105))) -3871833.4425726123)
  164.   (f2cl-lib:fset (f2cl-lib:fref c (105) ((1 105))) 18257.755474293175)
  165.   (f2cl-lib:fset (f2cl-lib:fref alfa (1) ((1 180))) -0.0044444444444444444)
  166.   (f2cl-lib:fset (f2cl-lib:fref alfa (2) ((1 180))) -9.220779220779222e-4)
  167.   (f2cl-lib:fset (f2cl-lib:fref alfa (3) ((1 180))) -8.848928848928851e-5)
  168.   (f2cl-lib:fset (f2cl-lib:fref alfa (4) ((1 180))) 1.6592768783244974e-4)
  169.   (f2cl-lib:fset (f2cl-lib:fref alfa (5) ((1 180))) 2.4669137274179292e-4)
  170.   (f2cl-lib:fset (f2cl-lib:fref alfa (6) ((1 180))) 2.6599558934625483e-4)
  171.   (f2cl-lib:fset (f2cl-lib:fref alfa (7) ((1 180))) 2.6182429706150095e-4)
  172.   (f2cl-lib:fset (f2cl-lib:fref alfa (8) ((1 180))) 2.4873043734465564e-4)
  173.   (f2cl-lib:fset (f2cl-lib:fref alfa (9) ((1 180))) 2.327210400832321e-4)
  174.   (f2cl-lib:fset (f2cl-lib:fref alfa (10) ((1 180))) 2.163624857123651e-4)
  175.   (f2cl-lib:fset (f2cl-lib:fref alfa (11) ((1 180))) 2.0073885876275238e-4)
  176.   (f2cl-lib:fset (f2cl-lib:fref alfa (12) ((1 180))) 1.8626763663754517e-4)
  177.   (f2cl-lib:fset (f2cl-lib:fref alfa (13) ((1 180))) 1.730607759178765e-4)
  178.   (f2cl-lib:fset (f2cl-lib:fref alfa (14) ((1 180))) 1.6109170592901578e-4)
  179.   (f2cl-lib:fset (f2cl-lib:fref alfa (15) ((1 180))) 1.5027477416090815e-4)
  180.   (f2cl-lib:fset (f2cl-lib:fref alfa (16) ((1 180))) 1.405034973912698e-4)
  181.   (f2cl-lib:fset (f2cl-lib:fref alfa (17) ((1 180))) 1.316688165459228e-4)
  182.   (f2cl-lib:fset (f2cl-lib:fref alfa (18) ((1 180))) 1.2366744559825327e-4)
  183.   (f2cl-lib:fset (f2cl-lib:fref alfa (19) ((1 180))) 1.1640527147473789e-4)
  184.   (f2cl-lib:fset (f2cl-lib:fref alfa (20) ((1 180))) 1.0979829837271338e-4)
  185.   (f2cl-lib:fset (f2cl-lib:fref alfa (21) ((1 180))) 1.0377241042299283e-4)
  186.   (f2cl-lib:fset (f2cl-lib:fref alfa (22) ((1 180))) 9.826260783693634e-5)
  187.   (f2cl-lib:fset (f2cl-lib:fref alfa (23) ((1 180))) 9.321205172495034e-5)
  188.   (f2cl-lib:fset (f2cl-lib:fref alfa (24) ((1 180))) 8.857108524787118e-5)
  189.   (f2cl-lib:fset (f2cl-lib:fref alfa (25) ((1 180))) 8.429631057157001e-5)
  190.   (f2cl-lib:fset (f2cl-lib:fref alfa (26) ((1 180))) 8.034975484077912e-5)
  191.   (f2cl-lib:fset (f2cl-lib:fref alfa (27) ((1 180))) 7.669813453592074e-5)
  192.   (f2cl-lib:fset (f2cl-lib:fref alfa (28) ((1 180))) 7.33122157481778e-5)
  193.   (f2cl-lib:fset (f2cl-lib:fref alfa (29) ((1 180))) 7.016626251631414e-5)
  194.   (f2cl-lib:fset (f2cl-lib:fref alfa (30) ((1 180))) 6.723756337901602e-5)
  195.   (f2cl-lib:fset (f2cl-lib:fref alfa (31) ((1 180))) 6.937355413545891e-4)
  196.   (f2cl-lib:fset (f2cl-lib:fref alfa (32) ((1 180))) 2.3224174518292165e-4)
  197.   (f2cl-lib:fset (f2cl-lib:fref alfa (33) ((1 180))) -1.4198627355669122e-5)
  198.   (f2cl-lib:fset (f2cl-lib:fref alfa (34) ((1 180))) -1.1644493167204864e-4)
  199.   (f2cl-lib:fset (f2cl-lib:fref alfa (35) ((1 180))) -1.508035580530488e-4)
  200.   (f2cl-lib:fset (f2cl-lib:fref alfa (36) ((1 180))) -1.5512192491809623e-4)
  201.   (f2cl-lib:fset (f2cl-lib:fref alfa (37) ((1 180))) -1.4680975664646556e-4)
  202.   (f2cl-lib:fset (f2cl-lib:fref alfa (38) ((1 180))) -1.3381550386749136e-4)
  203.   (f2cl-lib:fset (f2cl-lib:fref alfa (39) ((1 180))) -1.1974497568425405e-4)
  204.   (f2cl-lib:fset (f2cl-lib:fref alfa (40) ((1 180))) -1.0618431920797404e-4)
  205.   (f2cl-lib:fset (f2cl-lib:fref alfa (41) ((1 180))) -9.376995498911946e-5)
  206.   (f2cl-lib:fset (f2cl-lib:fref alfa (42) ((1 180))) -8.269230455881934e-5)
  207.   (f2cl-lib:fset (f2cl-lib:fref alfa (43) ((1 180))) -7.293743481552213e-5)
  208.   (f2cl-lib:fset (f2cl-lib:fref alfa (44) ((1 180))) -6.440423577210163e-5)
  209.   (f2cl-lib:fset (f2cl-lib:fref alfa (45) ((1 180))) -5.69611566009369e-5)
  210.   (f2cl-lib:fset (f2cl-lib:fref alfa (46) ((1 180))) -5.047310443035617e-5)
  211.   (f2cl-lib:fset (f2cl-lib:fref alfa (47) ((1 180))) -4.481348680088828e-5)
  212.   (f2cl-lib:fset (f2cl-lib:fref alfa (48) ((1 180))) -3.9868872771759895e-5)
  213.   (f2cl-lib:fset (f2cl-lib:fref alfa (49) ((1 180))) -3.554005329720425e-5)
  214.   (f2cl-lib:fset (f2cl-lib:fref alfa (50) ((1 180))) -3.174142566090225e-5)
  215.   (f2cl-lib:fset (f2cl-lib:fref alfa (51) ((1 180))) -2.8399679390417487e-5)
  216.   (f2cl-lib:fset (f2cl-lib:fref alfa (52) ((1 180))) -2.545227206348706e-5)
  217.   (f2cl-lib:fset (f2cl-lib:fref alfa (53) ((1 180))) -2.2845929716472457e-5)
  218.   (f2cl-lib:fset (f2cl-lib:fref alfa (54) ((1 180))) -2.0535275310648063e-5)
  219.   (f2cl-lib:fset (f2cl-lib:fref alfa (55) ((1 180))) -1.8481621762766612e-5)
  220.   (f2cl-lib:fset (f2cl-lib:fref alfa (56) ((1 180))) -1.665193300213938e-5)
  221.   (f2cl-lib:fset (f2cl-lib:fref alfa (57) ((1 180))) -1.501794129801195e-5)
  222.   (f2cl-lib:fset (f2cl-lib:fref alfa (58) ((1 180))) -1.3555403137904054e-5)
  223.   (f2cl-lib:fset (f2cl-lib:fref alfa (59) ((1 180))) -1.2243474647385814e-5)
  224.   (f2cl-lib:fset (f2cl-lib:fref alfa (60) ((1 180))) -1.1064188481130818e-5)
  225.   (f2cl-lib:fset (f2cl-lib:fref alfa (61) ((1 180))) -3.5421197145774386e-4)
  226.   (f2cl-lib:fset (f2cl-lib:fref alfa (62) ((1 180))) -1.5616126394515942e-4)
  227.   (f2cl-lib:fset (f2cl-lib:fref alfa (63) ((1 180))) 3.0446550359493646e-5)
  228.   (f2cl-lib:fset (f2cl-lib:fref alfa (64) ((1 180))) 1.3019865577324272e-4)
  229.   (f2cl-lib:fset (f2cl-lib:fref alfa (65) ((1 180))) 1.6747110669971227e-4)
  230.   (f2cl-lib:fset (f2cl-lib:fref alfa (66) ((1 180))) 1.7022258768359258e-4)
  231.   (f2cl-lib:fset (f2cl-lib:fref alfa (67) ((1 180))) 1.565014276085947e-4)
  232.   (f2cl-lib:fset (f2cl-lib:fref alfa (68) ((1 180))) 1.3633917097744513e-4)
  233.   (f2cl-lib:fset (f2cl-lib:fref alfa (69) ((1 180))) 1.1488669202982513e-4)
  234.   (f2cl-lib:fset (f2cl-lib:fref alfa (70) ((1 180))) 9.458690930346881e-5)
  235.   (f2cl-lib:fset (f2cl-lib:fref alfa (71) ((1 180))) 7.644984192508983e-5)
  236.   (f2cl-lib:fset (f2cl-lib:fref alfa (72) ((1 180))) 6.075703349651974e-5)
  237.   (f2cl-lib:fset (f2cl-lib:fref alfa (73) ((1 180))) 4.743942992905088e-5)
  238.   (f2cl-lib:fset (f2cl-lib:fref alfa (74) ((1 180))) 3.627575120053444e-5)
  239.   (f2cl-lib:fset (f2cl-lib:fref alfa (75) ((1 180))) 2.69939714979225e-5)
  240.   (f2cl-lib:fset (f2cl-lib:fref alfa (76) ((1 180))) 1.9321093824793925e-5)
  241.   (f2cl-lib:fset (f2cl-lib:fref alfa (77) ((1 180))) 1.300566747939632e-5)
  242.   (f2cl-lib:fset (f2cl-lib:fref alfa (78) ((1 180))) 7.826208667444966e-6)
  243.   (f2cl-lib:fset (f2cl-lib:fref alfa (79) ((1 180))) 3.592574858193515e-6)
  244.   (f2cl-lib:fset (f2cl-lib:fref alfa (80) ((1 180))) 1.4404004981425184e-7)
  245.   (f2cl-lib:fset (f2cl-lib:fref alfa (81) ((1 180))) -2.653967696979391e-6)
  246.   (f2cl-lib:fset (f2cl-lib:fref alfa (82) ((1 180))) -4.9134686709848585e-6)
  247.   (f2cl-lib:fset (f2cl-lib:fref alfa (83) ((1 180))) -6.727392960912483e-6)
  248.   (f2cl-lib:fset (f2cl-lib:fref alfa (84) ((1 180))) -8.172693796786579e-6)
  249.   (f2cl-lib:fset (f2cl-lib:fref alfa (85) ((1 180))) -9.313047150935612e-6)
  250.   (f2cl-lib:fset (f2cl-lib:fref alfa (86) ((1 180))) -1.0201141879801645e-5)
  251.   (f2cl-lib:fset (f2cl-lib:fref alfa (87) ((1 180))) -1.088059625105929e-5)
  252.   (f2cl-lib:fset (f2cl-lib:fref alfa (88) ((1 180))) -1.1387548150960355e-5)
  253.   (f2cl-lib:fset (f2cl-lib:fref alfa (89) ((1 180))) -1.1751967567455643e-5)
  254.   (f2cl-lib:fset (f2cl-lib:fref alfa (90) ((1 180))) -1.1998736487094415e-5)
  255.   (f2cl-lib:fset (f2cl-lib:fref alfa (91) ((1 180))) 3.781941992017729e-4)
  256.   (f2cl-lib:fset (f2cl-lib:fref alfa (92) ((1 180))) 2.0247195276181618e-4)
  257.   (f2cl-lib:fset (f2cl-lib:fref alfa (93) ((1 180))) -6.379385063188625e-5)
  258.   (f2cl-lib:fset (f2cl-lib:fref alfa (94) ((1 180))) -2.385982306030059e-4)
  259.   (f2cl-lib:fset (f2cl-lib:fref alfa (95) ((1 180))) -3.109162560273616e-4)
  260.   (f2cl-lib:fset (f2cl-lib:fref alfa (96) ((1 180))) -3.136801152475764e-4)
  261.   (f2cl-lib:fset (f2cl-lib:fref alfa (97) ((1 180))) -2.789502737913234e-4)
  262.   (f2cl-lib:fset (f2cl-lib:fref alfa (98) ((1 180))) -2.2856408261914143e-4)
  263.   (f2cl-lib:fset (f2cl-lib:fref alfa (99) ((1 180))) -1.7524528034084677e-4)
  264.   (f2cl-lib:fset (f2cl-lib:fref alfa (100) ((1 180))) -1.2554406306069035e-4)
  265.   (f2cl-lib:fset (f2cl-lib:fref alfa (101) ((1 180))) -8.229828728202083e-5)
  266.   (f2cl-lib:fset (f2cl-lib:fref alfa (102) ((1 180))) -4.628607305881165e-5)
  267.   (f2cl-lib:fset (f2cl-lib:fref alfa (103) ((1 180))) -1.7233430236696232e-5)
  268.   (f2cl-lib:fset (f2cl-lib:fref alfa (104) ((1 180))) 5.606904823046023e-6)
  269.   (f2cl-lib:fset (f2cl-lib:fref alfa (105) ((1 180))) 2.313954431482868e-5)
  270.   (f2cl-lib:fset (f2cl-lib:fref alfa (106) ((1 180))) 3.62642745856794e-5)
  271.   (f2cl-lib:fset (f2cl-lib:fref alfa (107) ((1 180))) 4.5800612449018885e-5)
  272.   (f2cl-lib:fset (f2cl-lib:fref alfa (108) ((1 180))) 5.245952949591141e-5)
  273.   (f2cl-lib:fset (f2cl-lib:fref alfa (109) ((1 180))) 5.683962085458154e-5)
  274.   (f2cl-lib:fset (f2cl-lib:fref alfa (110) ((1 180))) 5.943498203931041e-5)
  275.   (f2cl-lib:fset (f2cl-lib:fref alfa (111) ((1 180))) 6.064785275784218e-5)
  276.   (f2cl-lib:fset (f2cl-lib:fref alfa (112) ((1 180))) 6.080239077884365e-5)
  277.   (f2cl-lib:fset (f2cl-lib:fref alfa (113) ((1 180))) 6.015778945394604e-5)
  278.   (f2cl-lib:fset (f2cl-lib:fref alfa (114) ((1 180))) 5.891996573446985e-5)
  279.   (f2cl-lib:fset (f2cl-lib:fref alfa (115) ((1 180))) 5.725158237775931e-5)
  280.   (f2cl-lib:fset (f2cl-lib:fref alfa (116) ((1 180))) 5.528043755858526e-5)
  281.   (f2cl-lib:fset (f2cl-lib:fref alfa (117) ((1 180))) 5.310637738028802e-5)
  282.   (f2cl-lib:fset (f2cl-lib:fref alfa (118) ((1 180))) 5.080693020123257e-5)
  283.   (f2cl-lib:fset (f2cl-lib:fref alfa (119) ((1 180))) 4.844186476200949e-5)
  284.   (f2cl-lib:fset (f2cl-lib:fref alfa (120) ((1 180))) 4.6056858160747544e-5)
  285.   (f2cl-lib:fset (f2cl-lib:fref alfa (121) ((1 180))) -6.911413972882943e-4)
  286.   (f2cl-lib:fset (f2cl-lib:fref alfa (122) ((1 180))) -4.2997663305887196e-4)
  287.   (f2cl-lib:fset (f2cl-lib:fref alfa (123) ((1 180))) 1.8306773598003906e-4)
  288.   (f2cl-lib:fset (f2cl-lib:fref alfa (124) ((1 180))) 6.600881475420142e-4)
  289.   (f2cl-lib:fset (f2cl-lib:fref alfa (125) ((1 180))) 8.75964969951186e-4)
  290.   (f2cl-lib:fset (f2cl-lib:fref alfa (126) ((1 180))) 8.773352359582357e-4)
  291.   (f2cl-lib:fset (f2cl-lib:fref alfa (127) ((1 180))) 7.493695853789907e-4)
  292.   (f2cl-lib:fset (f2cl-lib:fref alfa (128) ((1 180))) 5.638323297569809e-4)
  293.   (f2cl-lib:fset (f2cl-lib:fref alfa (129) ((1 180))) 3.680593199714432e-4)
  294.   (f2cl-lib:fset (f2cl-lib:fref alfa (130) ((1 180))) 1.884645355144556e-4)
  295.   (f2cl-lib:fset (f2cl-lib:fref alfa (131) ((1 180))) 3.7066305766490415e-5)
  296.   (f2cl-lib:fset (f2cl-lib:fref alfa (132) ((1 180))) -8.285202202321371e-5)
  297.   (f2cl-lib:fset (f2cl-lib:fref alfa (133) ((1 180))) -1.72751952869173e-4)
  298.   (f2cl-lib:fset (f2cl-lib:fref alfa (134) ((1 180))) -2.3631487360587297e-4)
  299.   (f2cl-lib:fset (f2cl-lib:fref alfa (135) ((1 180))) -2.779661506949067e-4)
  300.   (f2cl-lib:fset (f2cl-lib:fref alfa (136) ((1 180))) -3.020795141554569e-4)
  301.   (f2cl-lib:fset (f2cl-lib:fref alfa (137) ((1 180))) -3.125947126438201e-4)
  302.   (f2cl-lib:fset (f2cl-lib:fref alfa (138) ((1 180))) -3.128725587580672e-4)
  303.   (f2cl-lib:fset (f2cl-lib:fref alfa (139) ((1 180))) -3.056780384663244e-4)
  304.   (f2cl-lib:fset (f2cl-lib:fref alfa (140) ((1 180))) -2.932264706145573e-4)
  305.   (f2cl-lib:fset (f2cl-lib:fref alfa (141) ((1 180))) -2.772556555829348e-4)
  306.   (f2cl-lib:fset (f2cl-lib:fref alfa (142) ((1 180))) -2.5910392846703173e-4)
  307.   (f2cl-lib:fset (f2cl-lib:fref alfa (143) ((1 180))) -2.3978401439648037e-4)
  308.   (f2cl-lib:fset (f2cl-lib:fref alfa (144) ((1 180))) -2.2004826004542286e-4)
  309.   (f2cl-lib:fset (f2cl-lib:fref alfa (145) ((1 180))) -2.0044391109497148e-4)
  310.   (f2cl-lib:fset (f2cl-lib:fref alfa (146) ((1 180))) -1.813586922109707e-4)
  311.   (f2cl-lib:fset (f2cl-lib:fref alfa (147) ((1 180))) -1.6305767447865746e-4)
  312.   (f2cl-lib:fset (f2cl-lib:fref alfa (148) ((1 180))) -1.4571267217520587e-4)
  313.   (f2cl-lib:fset (f2cl-lib:fref alfa (149) ((1 180))) -1.2942542198392462e-4)
  314.   (f2cl-lib:fset (f2cl-lib:fref alfa (150) ((1 180))) -1.1424569194244598e-4)
  315.   (f2cl-lib:fset (f2cl-lib:fref alfa (151) ((1 180))) 0.0019282196424877589)
  316.   (f2cl-lib:fset (f2cl-lib:fref alfa (152) ((1 180))) 0.0013559257630202223)
  317.   (f2cl-lib:fset (f2cl-lib:fref alfa (153) ((1 180))) -7.178580904213031e-4)
  318.   (f2cl-lib:fset (f2cl-lib:fref alfa (154) ((1 180))) -0.0025808480257527035)
  319.   (f2cl-lib:fset (f2cl-lib:fref alfa (155) ((1 180))) -0.003492711308261685)
  320.   (f2cl-lib:fset (f2cl-lib:fref alfa (156) ((1 180))) -0.003469862993409606)
  321.   (f2cl-lib:fset (f2cl-lib:fref alfa (157) ((1 180))) -0.002822852333513102)
  322.   (f2cl-lib:fset (f2cl-lib:fref alfa (158) ((1 180))) -0.0018810307640489134)
  323.   (f2cl-lib:fset (f2cl-lib:fref alfa (159) ((1 180))) -8.895317183839476e-4)
  324.   (f2cl-lib:fset (f2cl-lib:fref alfa (160) ((1 180))) 3.879121026310352e-6)
  325.   (f2cl-lib:fset (f2cl-lib:fref alfa (161) ((1 180))) 7.286885401196914e-4)
  326.   (f2cl-lib:fset (f2cl-lib:fref alfa (162) ((1 180))) 0.0012656637305345775)
  327.   (f2cl-lib:fset (f2cl-lib:fref alfa (163) ((1 180))) 0.0016251815837267443)
  328.   (f2cl-lib:fset (f2cl-lib:fref alfa (164) ((1 180))) 0.0018320315321637317)
  329.   (f2cl-lib:fset (f2cl-lib:fref alfa (165) ((1 180))) 0.001915883889905279)
  330.   (f2cl-lib:fset (f2cl-lib:fref alfa (166) ((1 180))) 0.0019058884675554615)
  331.   (f2cl-lib:fset (f2cl-lib:fref alfa (167) ((1 180))) 0.0018279898242182574)
  332.   (f2cl-lib:fset (f2cl-lib:fref alfa (168) ((1 180))) 0.0017038950642112155)
  333.   (f2cl-lib:fset (f2cl-lib:fref alfa (169) ((1 180))) 0.001550971271710977)
  334.   (f2cl-lib:fset (f2cl-lib:fref alfa (170) ((1 180))) 0.0013826142185227616)
  335.   (f2cl-lib:fset (f2cl-lib:fref alfa (171) ((1 180))) 0.0012088142423006478)
  336.   (f2cl-lib:fset (f2cl-lib:fref alfa (172) ((1 180))) 0.0010367653263834496)
  337.   (f2cl-lib:fset (f2cl-lib:fref alfa (173) ((1 180))) 8.714379180686191e-4)
  338.   (f2cl-lib:fset (f2cl-lib:fref alfa (174) ((1 180))) 7.16080155297701e-4)
  339.   (f2cl-lib:fset (f2cl-lib:fref alfa (175) ((1 180))) 5.726370025581293e-4)
  340.   (f2cl-lib:fset (f2cl-lib:fref alfa (176) ((1 180))) 4.420898194658023e-4)
  341.   (f2cl-lib:fset (f2cl-lib:fref alfa (177) ((1 180))) 3.247249485030906e-4)
  342.   (f2cl-lib:fset (f2cl-lib:fref alfa (178) ((1 180))) 2.2034204273024657e-4)
  343.   (f2cl-lib:fset (f2cl-lib:fref alfa (179) ((1 180))) 1.284128984013539e-4)
  344.   (f2cl-lib:fset (f2cl-lib:fref alfa (180) ((1 180))) 4.8200592455209545e-5)
  345.   (f2cl-lib:fset (f2cl-lib:fref beta (1) ((1 210))) 0.017998872141355333)
  346.   (f2cl-lib:fset (f2cl-lib:fref beta (2) ((1 210))) 0.005599649110643881)
  347.   (f2cl-lib:fset (f2cl-lib:fref beta (3) ((1 210))) 0.002885014022311328)
  348.   (f2cl-lib:fset (f2cl-lib:fref beta (4) ((1 210))) 0.0018009660676105395)
  349.   (f2cl-lib:fset (f2cl-lib:fref beta (5) ((1 210))) 0.001247531105891992)
  350.   (f2cl-lib:fset (f2cl-lib:fref beta (6) ((1 210))) 9.228788765729384e-4)
  351.   (f2cl-lib:fset (f2cl-lib:fref beta (7) ((1 210))) 7.144304217272874e-4)
  352.   (f2cl-lib:fset (f2cl-lib:fref beta (8) ((1 210))) 5.71787281789705e-4)
  353.   (f2cl-lib:fset (f2cl-lib:fref beta (9) ((1 210))) 4.6943100760648154e-4)
  354.   (f2cl-lib:fset (f2cl-lib:fref beta (10) ((1 210))) 3.932328354629167e-4)
  355.   (f2cl-lib:fset (f2cl-lib:fref beta (11) ((1 210))) 3.348188893182977e-4)
  356.   (f2cl-lib:fset (f2cl-lib:fref beta (12) ((1 210))) 2.8895214849575157e-4)
  357.   (f2cl-lib:fset (f2cl-lib:fref beta (13) ((1 210))) 2.5221161554957333e-4)
  358.   (f2cl-lib:fset (f2cl-lib:fref beta (14) ((1 210))) 2.2228058079888333e-4)
  359.   (f2cl-lib:fset (f2cl-lib:fref beta (15) ((1 210))) 1.9754183803306255e-4)
  360.   (f2cl-lib:fset (f2cl-lib:fref beta (16) ((1 210))) 1.7683685501971802e-4)
  361.   (f2cl-lib:fset (f2cl-lib:fref beta (17) ((1 210))) 1.593168996618211e-4)
  362.   (f2cl-lib:fset (f2cl-lib:fref beta (18) ((1 210))) 1.44347930197334e-4)
  363.   (f2cl-lib:fset (f2cl-lib:fref beta (19) ((1 210))) 1.314480681199654e-4)
  364.   (f2cl-lib:fset (f2cl-lib:fref beta (20) ((1 210))) 1.2024544494930292e-4)
  365.   (f2cl-lib:fset (f2cl-lib:fref beta (21) ((1 210))) 1.1044914450459939e-4)
  366.   (f2cl-lib:fset (f2cl-lib:fref beta (22) ((1 210))) 1.0182877074056726e-4)
  367.   (f2cl-lib:fset (f2cl-lib:fref beta (23) ((1 210))) 9.419982242042375e-5)
  368.   (f2cl-lib:fset (f2cl-lib:fref beta (24) ((1 210))) 8.741305457538344e-5)
  369.   (f2cl-lib:fset (f2cl-lib:fref beta (25) ((1 210))) 8.134662621628016e-5)
  370.   (f2cl-lib:fset (f2cl-lib:fref beta (26) ((1 210))) 7.590022696462193e-5)
  371.   (f2cl-lib:fset (f2cl-lib:fref beta (27) ((1 210))) 7.099063006341535e-5)
  372.   (f2cl-lib:fset (f2cl-lib:fref beta (28) ((1 210))) 6.654828748424683e-5)
  373.   (f2cl-lib:fset (f2cl-lib:fref beta (29) ((1 210))) 6.251469589692751e-5)
  374.   (f2cl-lib:fset (f2cl-lib:fref beta (30) ((1 210))) 5.884033944262518e-5)
  375.   (f2cl-lib:fset (f2cl-lib:fref beta (31) ((1 210))) -0.0014928295321342917)
  376.   (f2cl-lib:fset (f2cl-lib:fref beta (32) ((1 210))) -8.782047095463893e-4)
  377.   (f2cl-lib:fset (f2cl-lib:fref beta (33) ((1 210))) -5.029165495720346e-4)
  378.   (f2cl-lib:fset (f2cl-lib:fref beta (34) ((1 210))) -2.948221385127461e-4)
  379.   (f2cl-lib:fset (f2cl-lib:fref beta (35) ((1 210))) -1.7546399697078285e-4)
  380.   (f2cl-lib:fset (f2cl-lib:fref beta (36) ((1 210))) -1.0400855046081643e-4)
  381.   (f2cl-lib:fset (f2cl-lib:fref beta (37) ((1 210))) -5.96141953046458e-5)
  382.   (f2cl-lib:fset (f2cl-lib:fref beta (38) ((1 210))) -3.120389290760983e-5)
  383.   (f2cl-lib:fset (f2cl-lib:fref beta (39) ((1 210))) -1.2608973598023004e-5)
  384.   (f2cl-lib:fset (f2cl-lib:fref beta (40) ((1 210))) -2.4289260857573036e-7)
  385.   (f2cl-lib:fset (f2cl-lib:fref beta (41) ((1 210))) 8.059961654142735e-6)
  386.   (f2cl-lib:fset (f2cl-lib:fref beta (42) ((1 210))) 1.3650700926214737e-5)
  387.   (f2cl-lib:fset (f2cl-lib:fref beta (43) ((1 210))) 1.7396412547292628e-5)
  388.   (f2cl-lib:fset (f2cl-lib:fref beta (44) ((1 210))) 1.9867297884213378e-5)
  389.   (f2cl-lib:fset (f2cl-lib:fref beta (45) ((1 210))) 2.1446326379082264e-5)
  390.   (f2cl-lib:fset (f2cl-lib:fref beta (46) ((1 210))) 2.2395465923245655e-5)
  391.   (f2cl-lib:fset (f2cl-lib:fref beta (47) ((1 210))) 2.289677838147126e-5)
  392.   (f2cl-lib:fset (f2cl-lib:fref beta (48) ((1 210))) 2.307853898111778e-5)
  393.   (f2cl-lib:fset (f2cl-lib:fref beta (49) ((1 210))) 2.3032197608090915e-5)
  394.   (f2cl-lib:fset (f2cl-lib:fref beta (50) ((1 210))) 2.2823607372034874e-5)
  395.   (f2cl-lib:fset (f2cl-lib:fref beta (51) ((1 210))) 2.250058811052924e-5)
  396.   (f2cl-lib:fset (f2cl-lib:fref beta (52) ((1 210))) 2.2098101536199147e-5)
  397.   (f2cl-lib:fset (f2cl-lib:fref beta (53) ((1 210))) 2.1641842744810388e-5)
  398.   (f2cl-lib:fset (f2cl-lib:fref beta (54) ((1 210))) 2.1150764925622084e-5)
  399.   (f2cl-lib:fset (f2cl-lib:fref beta (55) ((1 210))) 2.0638874978217076e-5)
  400.   (f2cl-lib:fset (f2cl-lib:fref beta (56) ((1 210))) 2.011652419970817e-5)
  401.   (f2cl-lib:fset (f2cl-lib:fref beta (57) ((1 210))) 1.9591345014117925e-5)
  402.   (f2cl-lib:fset (f2cl-lib:fref beta (58) ((1 210))) 1.9068936791043675e-5)
  403.   (f2cl-lib:fset (f2cl-lib:fref beta (59) ((1 210))) 1.855337196416367e-5)
  404.   (f2cl-lib:fset (f2cl-lib:fref beta (60) ((1 210))) 1.8047572225967418e-5)
  405.   (f2cl-lib:fset (f2cl-lib:fref beta (61) ((1 210))) 5.522130767212928e-4)
  406.   (f2cl-lib:fset (f2cl-lib:fref beta (62) ((1 210))) 4.4793258155238475e-4)
  407.   (f2cl-lib:fset (f2cl-lib:fref beta (63) ((1 210))) 2.795206539920206e-4)
  408.   (f2cl-lib:fset (f2cl-lib:fref beta (64) ((1 210))) 1.5246815619844659e-4)
  409.   (f2cl-lib:fset (f2cl-lib:fref beta (65) ((1 210))) 6.932711056570436e-5)
  410.   (f2cl-lib:fset (f2cl-lib:fref beta (66) ((1 210))) 1.7625868306999143e-5)
  411.   (f2cl-lib:fset (f2cl-lib:fref beta (67) ((1 210))) -1.3574499634326914e-5)
  412.   (f2cl-lib:fset (f2cl-lib:fref beta (68) ((1 210))) -3.1797241335042714e-5)
  413.   (f2cl-lib:fset (f2cl-lib:fref beta (69) ((1 210))) -4.188618616966934e-5)
  414.   (f2cl-lib:fset (f2cl-lib:fref beta (70) ((1 210))) -4.690048893791411e-5)
  415.   (f2cl-lib:fset (f2cl-lib:fref beta (71) ((1 210))) -4.8766544741378737e-5)
  416.   (f2cl-lib:fset (f2cl-lib:fref beta (72) ((1 210))) -4.8701003118673514e-5)
  417.   (f2cl-lib:fset (f2cl-lib:fref beta (73) ((1 210))) -4.747556208900867e-5)
  418.   (f2cl-lib:fset (f2cl-lib:fref beta (74) ((1 210))) -4.558130581386285e-5)
  419.   (f2cl-lib:fset (f2cl-lib:fref beta (75) ((1 210))) -4.3330964451126613e-5)
  420.   (f2cl-lib:fset (f2cl-lib:fref beta (76) ((1 210))) -4.0923019315775044e-5)
  421.   (f2cl-lib:fset (f2cl-lib:fref beta (77) ((1 210))) -3.848226386032213e-5)
  422.   (f2cl-lib:fset (f2cl-lib:fref beta (78) ((1 210))) -3.6085716753541053e-5)
  423.   (f2cl-lib:fset (f2cl-lib:fref beta (79) ((1 210))) -3.3779330612336744e-5)
  424.   (f2cl-lib:fset (f2cl-lib:fref beta (80) ((1 210))) -3.158885607721097e-5)
  425.   (f2cl-lib:fset (f2cl-lib:fref beta (81) ((1 210))) -2.952695617508073e-5)
  426.   (f2cl-lib:fset (f2cl-lib:fref beta (82) ((1 210))) -2.7597891482833575e-5)
  427.   (f2cl-lib:fset (f2cl-lib:fref beta (83) ((1 210))) -2.580061746668837e-5)
  428.   (f2cl-lib:fset (f2cl-lib:fref beta (84) ((1 210))) -2.413083567612802e-5)
  429.   (f2cl-lib:fset (f2cl-lib:fref beta (85) ((1 210))) -2.2582350951834607e-5)
  430.   (f2cl-lib:fset (f2cl-lib:fref beta (86) ((1 210))) -2.1147965676891295e-5)
  431.   (f2cl-lib:fset (f2cl-lib:fref beta (87) ((1 210))) -1.982006388852949e-5)
  432.   (f2cl-lib:fset (f2cl-lib:fref beta (88) ((1 210))) -1.8590987080106508e-5)
  433.   (f2cl-lib:fset (f2cl-lib:fref beta (89) ((1 210))) -1.7453269984421022e-5)
  434.   (f2cl-lib:fset (f2cl-lib:fref beta (90) ((1 210))) -1.6399782385449801e-5)
  435.   (f2cl-lib:fset (f2cl-lib:fref beta (91) ((1 210))) -4.746177965599597e-4)
  436.   (f2cl-lib:fset (f2cl-lib:fref beta (92) ((1 210))) -4.778645671473215e-4)
  437.   (f2cl-lib:fset (f2cl-lib:fref beta (93) ((1 210))) -3.2039022806703765e-4)
  438.   (f2cl-lib:fset (f2cl-lib:fref beta (94) ((1 210))) -1.6110501611996228e-4)
  439.   (f2cl-lib:fset (f2cl-lib:fref beta (95) ((1 210))) -4.2577810128543525e-5)
  440.   (f2cl-lib:fset (f2cl-lib:fref beta (96) ((1 210))) 3.445712942949675e-5)
  441.   (f2cl-lib:fset (f2cl-lib:fref beta (97) ((1 210))) 7.97092684075675e-5)
  442.   (f2cl-lib:fset (f2cl-lib:fref beta (98) ((1 210))) 1.031382367082722e-4)
  443.   (f2cl-lib:fset (f2cl-lib:fref beta (99) ((1 210))) 1.1246677526220417e-4)
  444.   (f2cl-lib:fset (f2cl-lib:fref beta (100) ((1 210))) 1.1310364210848138e-4)
  445.   (f2cl-lib:fset (f2cl-lib:fref beta (101) ((1 210))) 1.0865163484877427e-4)
  446.   (f2cl-lib:fset (f2cl-lib:fref beta (102) ((1 210))) 1.0143795159766197e-4)
  447.   (f2cl-lib:fset (f2cl-lib:fref beta (103) ((1 210))) 9.292983965933639e-5)
  448.   (f2cl-lib:fset (f2cl-lib:fref beta (104) ((1 210))) 8.402931330160901e-5)
  449.   (f2cl-lib:fset (f2cl-lib:fref beta (105) ((1 210))) 7.527279913491342e-5)
  450.   (f2cl-lib:fset (f2cl-lib:fref beta (106) ((1 210))) 6.696325219757311e-5)
  451.   (f2cl-lib:fset (f2cl-lib:fref beta (107) ((1 210))) 5.925645473231947e-5)
  452.   (f2cl-lib:fset (f2cl-lib:fref beta (108) ((1 210))) 5.221693088269757e-5)
  453.   (f2cl-lib:fset (f2cl-lib:fref beta (109) ((1 210))) 4.585394851653607e-5)
  454.   (f2cl-lib:fset (f2cl-lib:fref beta (110) ((1 210))) 4.014455138914868e-5)
  455.   (f2cl-lib:fset (f2cl-lib:fref beta (111) ((1 210))) 3.5048173003132804e-5)
  456.   (f2cl-lib:fset (f2cl-lib:fref beta (112) ((1 210))) 3.0515799503434665e-5)
  457.   (f2cl-lib:fset (f2cl-lib:fref beta (113) ((1 210))) 2.6495611995051604e-5)
  458.   (f2cl-lib:fset (f2cl-lib:fref beta (114) ((1 210))) 2.293636336909982e-5)
  459.   (f2cl-lib:fset (f2cl-lib:fref beta (115) ((1 210))) 1.9789305666402165e-5)
  460.   (f2cl-lib:fset (f2cl-lib:fref beta (116) ((1 210))) 1.7009198463641262e-5)
  461.   (f2cl-lib:fset (f2cl-lib:fref beta (117) ((1 210))) 1.4554742826152403e-5)
  462.   (f2cl-lib:fset (f2cl-lib:fref beta (118) ((1 210))) 1.2388664099587843e-5)
  463.   (f2cl-lib:fset (f2cl-lib:fref beta (119) ((1 210))) 1.0477587607658323e-5)
  464.   (f2cl-lib:fset (f2cl-lib:fref beta (120) ((1 210))) 8.791799549784793e-6)
  465.   (f2cl-lib:fset (f2cl-lib:fref beta (121) ((1 210))) 7.364658105725783e-4)
  466.   (f2cl-lib:fset (f2cl-lib:fref beta (122) ((1 210))) 8.72790805146194e-4)
  467.   (f2cl-lib:fset (f2cl-lib:fref beta (123) ((1 210))) 6.22614862573135e-4)
  468.   (f2cl-lib:fset (f2cl-lib:fref beta (124) ((1 210))) 2.859981541943042e-4)
  469.   (f2cl-lib:fset (f2cl-lib:fref beta (125) ((1 210))) 3.8473767287936605e-6)
  470.   (f2cl-lib:fset (f2cl-lib:fref beta (126) ((1 210))) -1.8790600363697157e-4)
  471.   (f2cl-lib:fset (f2cl-lib:fref beta (127) ((1 210))) -2.9760364659455457e-4)
  472.   (f2cl-lib:fset (f2cl-lib:fref beta (128) ((1 210))) -3.459981268326564e-4)
  473.   (f2cl-lib:fset (f2cl-lib:fref beta (129) ((1 210))) -3.533824709160377e-4)
  474.   (f2cl-lib:fset (f2cl-lib:fref beta (130) ((1 210))) -3.357156357750488e-4)
  475.   (f2cl-lib:fset (f2cl-lib:fref beta (131) ((1 210))) -3.043211247890398e-4)
  476.   (f2cl-lib:fset (f2cl-lib:fref beta (132) ((1 210))) -2.667227230476128e-4)
  477.   (f2cl-lib:fset (f2cl-lib:fref beta (133) ((1 210))) -2.2765421412281958e-4)
  478.   (f2cl-lib:fset (f2cl-lib:fref beta (134) ((1 210))) -1.8992261185456238e-4)
  479.   (f2cl-lib:fset (f2cl-lib:fref beta (135) ((1 210))) -1.5505891859909388e-4)
  480.   (f2cl-lib:fset (f2cl-lib:fref beta (136) ((1 210))) -1.2377824076187363e-4)
  481.   (f2cl-lib:fset (f2cl-lib:fref beta (137) ((1 210))) -9.629261477176442e-5)
  482.   (f2cl-lib:fset (f2cl-lib:fref beta (138) ((1 210))) -7.251783277144254e-5)
  483.   (f2cl-lib:fset (f2cl-lib:fref beta (139) ((1 210))) -5.220700288956338e-5)
  484.   (f2cl-lib:fset (f2cl-lib:fref beta (140) ((1 210))) -3.5034775051190054e-5)
  485.   (f2cl-lib:fset (f2cl-lib:fref beta (141) ((1 210))) -2.0648976103555175e-5)
  486.   (f2cl-lib:fset (f2cl-lib:fref beta (142) ((1 210))) -8.70106096849767e-6)
  487.   (f2cl-lib:fset (f2cl-lib:fref beta (143) ((1 210))) 1.1369868667510027e-6)
  488.   (f2cl-lib:fset (f2cl-lib:fref beta (144) ((1 210))) 9.164264741227788e-6)
  489.   (f2cl-lib:fset (f2cl-lib:fref beta (145) ((1 210))) 1.5647778542887264e-5)
  490.   (f2cl-lib:fset (f2cl-lib:fref beta (146) ((1 210))) 2.0822362948246687e-5)
  491.   (f2cl-lib:fset (f2cl-lib:fref beta (147) ((1 210))) 2.4892338100459516e-5)
  492.   (f2cl-lib:fset (f2cl-lib:fref beta (148) ((1 210))) 2.8034050957414636e-5)
  493.   (f2cl-lib:fset (f2cl-lib:fref beta (149) ((1 210))) 3.0398777462986193e-5)
  494.   (f2cl-lib:fset (f2cl-lib:fref beta (150) ((1 210))) 3.2115673140670065e-5)
  495.   (f2cl-lib:fset (f2cl-lib:fref beta (151) ((1 210))) -0.0018018219196388571)
  496.   (f2cl-lib:fset (f2cl-lib:fref beta (152) ((1 210))) -0.0024340296293804253)
  497.   (f2cl-lib:fset (f2cl-lib:fref beta (153) ((1 210))) -0.0018342266354985683)
  498.   (f2cl-lib:fset (f2cl-lib:fref beta (154) ((1 210))) -7.622045963540097e-4)
  499.   (f2cl-lib:fset (f2cl-lib:fref beta (155) ((1 210))) 2.3907947525692724e-4)
  500.   (f2cl-lib:fset (f2cl-lib:fref beta (156) ((1 210))) 9.49266117176881e-4)
  501.   (f2cl-lib:fset (f2cl-lib:fref beta (157) ((1 210))) 0.0013446744970154036)
  502.   (f2cl-lib:fset (f2cl-lib:fref beta (158) ((1 210))) 0.0014845749525944918)
  503.   (f2cl-lib:fset (f2cl-lib:fref beta (159) ((1 210))) 0.001447323398306176)
  504.   (f2cl-lib:fset (f2cl-lib:fref beta (160) ((1 210))) 0.0013026826128565718)
  505.   (f2cl-lib:fset (f2cl-lib:fref beta (161) ((1 210))) 0.0011035159737564268)
  506.   (f2cl-lib:fset (f2cl-lib:fref beta (162) ((1 210))) 8.860474404197918e-4)
  507.   (f2cl-lib:fset (f2cl-lib:fref beta (163) ((1 210))) 6.730732081656655e-4)
  508.   (f2cl-lib:fset (f2cl-lib:fref beta (164) ((1 210))) 4.776038728565824e-4)
  509.   (f2cl-lib:fset (f2cl-lib:fref beta (165) ((1 210))) 3.0599192635878936e-4)
  510.   (f2cl-lib:fset (f2cl-lib:fref beta (166) ((1 210))) 1.6031569459472167e-4)
  511.   (f2cl-lib:fset (f2cl-lib:fref beta (167) ((1 210))) 4.0074955527061334e-5)
  512.   (f2cl-lib:fset (f2cl-lib:fref beta (168) ((1 210))) -5.666074616352517e-5)
  513.   (f2cl-lib:fset (f2cl-lib:fref beta (169) ((1 210))) -1.3250618677298265e-4)
  514.   (f2cl-lib:fset (f2cl-lib:fref beta (170) ((1 210))) -1.9029618798961406e-4)
  515.   (f2cl-lib:fset (f2cl-lib:fref beta (171) ((1 210))) -2.328114503769374e-4)
  516.   (f2cl-lib:fset (f2cl-lib:fref beta (172) ((1 210))) -2.626288114646688e-4)
  517.   (f2cl-lib:fset (f2cl-lib:fref beta (173) ((1 210))) -2.8205046986759863e-4)
  518.   (f2cl-lib:fset (f2cl-lib:fref beta (174) ((1 210))) -2.9308156319286116e-4)
  519.   (f2cl-lib:fset (f2cl-lib:fref beta (175) ((1 210))) -2.974359621763166e-4)
  520.   (f2cl-lib:fset (f2cl-lib:fref beta (176) ((1 210))) -2.9655733423934805e-4)
  521.   (f2cl-lib:fset (f2cl-lib:fref beta (177) ((1 210))) -2.9164736331209085e-4)
  522.   (f2cl-lib:fset (f2cl-lib:fref beta (178) ((1 210))) -2.836962038377342e-4)
  523.   (f2cl-lib:fset (f2cl-lib:fref beta (179) ((1 210))) -2.7351231709567336e-4)
  524.   (f2cl-lib:fset (f2cl-lib:fref beta (180) ((1 210))) -2.617501558067686e-4)
  525.   (f2cl-lib:fset (f2cl-lib:fref beta (181) ((1 210))) 0.006385858912120509)
  526.   (f2cl-lib:fset (f2cl-lib:fref beta (182) ((1 210))) 0.009623742158063778)
  527.   (f2cl-lib:fset (f2cl-lib:fref beta (183) ((1 210))) 0.007618780612070011)
  528.   (f2cl-lib:fset (f2cl-lib:fref beta (184) ((1 210))) 0.002832190555456281)
  529.   (f2cl-lib:fset (f2cl-lib:fref beta (185) ((1 210))) -0.002098413520127201)
  530.   (f2cl-lib:fset (f2cl-lib:fref beta (186) ((1 210))) -0.005738267642166265)
  531.   (f2cl-lib:fset (f2cl-lib:fref beta (187) ((1 210))) -0.0077080424449541465)
  532.   (f2cl-lib:fset (f2cl-lib:fref beta (188) ((1 210))) -0.008210116922648444)
  533.   (f2cl-lib:fset (f2cl-lib:fref beta (189) ((1 210))) -0.007658245203469054)
  534.   (f2cl-lib:fset (f2cl-lib:fref beta (190) ((1 210))) -0.006472097293910452)
  535.   (f2cl-lib:fset (f2cl-lib:fref beta (191) ((1 210))) -0.004991324120049665)
  536.   (f2cl-lib:fset (f2cl-lib:fref beta (192) ((1 210))) -0.003456122897131333)
  537.   (f2cl-lib:fset (f2cl-lib:fref beta (193) ((1 210))) -0.002017855800141708)
  538.   (f2cl-lib:fset (f2cl-lib:fref beta (194) ((1 210))) -7.594306867819615e-4)
  539.   (f2cl-lib:fset (f2cl-lib:fref beta (195) ((1 210))) 2.8417363152385916e-4)
  540.   (f2cl-lib:fset (f2cl-lib:fref beta (196) ((1 210))) 0.0011089166758633739)
  541.   (f2cl-lib:fset (f2cl-lib:fref beta (197) ((1 210))) 0.0017290149387272878)
  542.   (f2cl-lib:fset (f2cl-lib:fref beta (198) ((1 210))) 0.0021681259080268468)
  543.   (f2cl-lib:fset (f2cl-lib:fref beta (199) ((1 210))) 0.0024535771049453976)
  544.   (f2cl-lib:fset (f2cl-lib:fref beta (200) ((1 210))) 0.0026128182105833483)
  545.   (f2cl-lib:fset (f2cl-lib:fref beta (201) ((1 210))) 0.002671410396562769)
  546.   (f2cl-lib:fset (f2cl-lib:fref beta (202) ((1 210))) 0.0026520307339598045)
  547.   (f2cl-lib:fset (f2cl-lib:fref beta (203) ((1 210))) 0.002574116528772873)
  548.   (f2cl-lib:fset (f2cl-lib:fref beta (204) ((1 210))) 0.0024538912623609443)
  549.   (f2cl-lib:fset (f2cl-lib:fref beta (205) ((1 210))) 0.002304600580717955)
  550.   (f2cl-lib:fset (f2cl-lib:fref beta (206) ((1 210))) 0.0021368483768671267)
  551.   (f2cl-lib:fset (f2cl-lib:fref beta (207) ((1 210))) 0.001958965284788709)
  552.   (f2cl-lib:fset (f2cl-lib:fref beta (208) ((1 210))) 0.0017773700867945441)
  553.   (f2cl-lib:fset (f2cl-lib:fref beta (209) ((1 210))) 0.0015969028076583908)
  554.   (f2cl-lib:fset (f2cl-lib:fref beta (210) ((1 210))) 0.0014211197566443854)
  555.   (f2cl-lib:fset (f2cl-lib:fref gama (1) ((1 30))) 0.6299605249474367)
  556.   (f2cl-lib:fset (f2cl-lib:fref gama (2) ((1 30))) 0.25198420997897464)
  557.   (f2cl-lib:fset (f2cl-lib:fref gama (3) ((1 30))) 0.15479030041565586)
  558.   (f2cl-lib:fset (f2cl-lib:fref gama (4) ((1 30))) 0.11071306241615902)
  559.   (f2cl-lib:fset (f2cl-lib:fref gama (5) ((1 30))) 0.08573093955273949)
  560.   (f2cl-lib:fset (f2cl-lib:fref gama (6) ((1 30))) 0.06971613169586843)
  561.   (f2cl-lib:fset (f2cl-lib:fref gama (7) ((1 30))) 0.05860856718937136)
  562.   (f2cl-lib:fset (f2cl-lib:fref gama (8) ((1 30))) 0.050469887353631074)
  563.   (f2cl-lib:fset (f2cl-lib:fref gama (9) ((1 30))) 0.04426005806891548)
  564.   (f2cl-lib:fset (f2cl-lib:fref gama (10) ((1 30))) 0.039372066154350994)
  565.   (f2cl-lib:fset (f2cl-lib:fref gama (11) ((1 30))) 0.03542831959244554)
  566.   (f2cl-lib:fset (f2cl-lib:fref gama (12) ((1 30))) 0.032181885750209825)
  567.   (f2cl-lib:fset (f2cl-lib:fref gama (13) ((1 30))) 0.029464624079115768)
  568.   (f2cl-lib:fset (f2cl-lib:fref gama (14) ((1 30))) 0.027158167711293448)
  569.   (f2cl-lib:fset (f2cl-lib:fref gama (15) ((1 30))) 0.025176827297386177)
  570.   (f2cl-lib:fset (f2cl-lib:fref gama (16) ((1 30))) 0.02345707553060789)
  571.   (f2cl-lib:fset (f2cl-lib:fref gama (17) ((1 30))) 0.021950839013490722)
  572.   (f2cl-lib:fset (f2cl-lib:fref gama (18) ((1 30))) 0.020621082823564625)
  573.   (f2cl-lib:fset (f2cl-lib:fref gama (19) ((1 30))) 0.019438824089788084)
  574.   (f2cl-lib:fset (f2cl-lib:fref gama (20) ((1 30))) 0.018381063380068317)
  575.   (f2cl-lib:fset (f2cl-lib:fref gama (21) ((1 30))) 0.017429321323196318)
  576.   (f2cl-lib:fset (f2cl-lib:fref gama (22) ((1 30))) 0.016568583778661234)
  577.   (f2cl-lib:fset (f2cl-lib:fref gama (23) ((1 30))) 0.015786528598791844)
  578.   (f2cl-lib:fset (f2cl-lib:fref gama (24) ((1 30))) 0.01507295014940956)
  579.   (f2cl-lib:fset (f2cl-lib:fref gama (25) ((1 30))) 0.014419325083995465)
  580.   (f2cl-lib:fset (f2cl-lib:fref gama (26) ((1 30))) 0.013818480573534178)
  581.   (f2cl-lib:fset (f2cl-lib:fref gama (27) ((1 30))) 0.013264337899427657)
  582.   (f2cl-lib:fset (f2cl-lib:fref gama (28) ((1 30))) 0.012751712197049864)
  583.   (f2cl-lib:fset (f2cl-lib:fref gama (29) ((1 30))) 0.012276154531876278)
  584.   (f2cl-lib:fset (f2cl-lib:fref gama (30) ((1 30))) 0.01183382623984824)
  585.   (defun zunhj
  586.          (zr zi fnu ipmtr tol phir phii argr argi zeta1r zeta1i zeta2r zeta2i
  587.           asumr asumi bsumr bsumi)
  588.     (declare (type f2cl-lib:integer4 ipmtr)
  589.              (type double-float zr zi fnu tol phir phii argr argi zeta1r zeta1i
  590.               zeta2r zeta2i asumr asumi bsumr bsumi))
  591.     (prog ((ap (make-array 30 :element-type 'double-float))
  592.            (pr (make-array 30 :element-type 'double-float))
  593.            (pi_ (make-array 30 :element-type 'double-float))
  594.            (upr (make-array 14 :element-type 'double-float))
  595.            (upi (make-array 14 :element-type 'double-float))
  596.            (crr (make-array 14 :element-type 'double-float))
  597.            (cri (make-array 14 :element-type 'double-float))
  598.            (drr (make-array 14 :element-type 'double-float))
  599.            (dri (make-array 14 :element-type 'double-float)) (ias 0) (ibs 0)
  600.            (is 0) (j 0) (jr 0) (ju 0) (k 0) (kmax 0) (kp1 0) (ks 0) (l 0)
  601.            (lr 0) (lrp1 0) (l1 0) (l2 0) (m 0) (idum 0) (ang 0.0) (atol 0.0)
  602.            (aw2 0.0) (azth 0.0) (btol 0.0) (fn13 0.0) (fn23 0.0) (pp 0.0)
  603.            (przthi 0.0) (przthr 0.0) (ptfni 0.0) (ptfnr 0.0) (raw 0.0)
  604.            (raw2 0.0) (razth 0.0) (rfnu 0.0) (rfnu2 0.0) (rfn13 0.0)
  605.            (rtzti 0.0) (rtztr 0.0) (rzthi 0.0) (rzthr 0.0) (sti 0.0) (str 0.0)
  606.            (sumai 0.0) (sumar 0.0) (sumbi 0.0) (sumbr 0.0) (test 0.0)
  607.            (tfni 0.0) (tfnr 0.0) (tzai 0.0) (tzar 0.0) (t2i 0.0) (t2r 0.0)
  608.            (wi 0.0) (wr 0.0) (w2i 0.0) (w2r 0.0) (zai 0.0) (zar 0.0) (zbi 0.0)
  609.            (zbr 0.0) (zci 0.0) (zcr 0.0) (zetai 0.0) (zetar 0.0) (zthi 0.0)
  610.            (zthr 0.0) (ac 0.0))
  611.       (declare (type (simple-array double-float (14)) upr upi drr dri crr cri)
  612.                (type (simple-array double-float (30)) pr pi_ ap)
  613.                (type double-float ac zthr zthi zetar zetai zcr zci zbr zbi zar
  614.                 zai w2r w2i wr wi t2r t2i tzar tzai tfnr tfni test sumbr sumbi
  615.                 sumar sumai str sti rzthr rzthi rtztr rtzti rfn13 rfnu2 rfnu
  616.                 razth raw2 raw ptfnr ptfni przthr przthi pp fn23 fn13 btol azth
  617.                 aw2 atol ang)
  618.                (type f2cl-lib:integer4 idum m l2 l1 lrp1 lr l ks kp1 kmax k ju
  619.                 jr j is ibs ias))
  620.       (setf rfnu (/ 1.0 fnu))
  621.       (setf test (* (f2cl-lib:d1mach 1) 1000.0))
  622.       (setf ac (* fnu test))
  623.       (if (or (> (abs zr) ac) (> (abs zi) ac)) (go label15))
  624.       (setf zeta1r (+ (* 2.0 (abs (f2cl-lib:flog test))) fnu))
  625.       (setf zeta1i 0.0)
  626.       (setf zeta2r fnu)
  627.       (setf zeta2i 0.0)
  628.       (setf phir 1.0)
  629.       (setf phii 0.0)
  630.       (setf argr 1.0)
  631.       (setf argi 0.0)
  632.       (go end_label)
  633.      label15
  634.       (setf zbr (* zr rfnu))
  635.       (setf zbi (* zi rfnu))
  636.       (setf rfnu2 (* rfnu rfnu))
  637.       (setf fn13 (expt fnu ex1))
  638.       (setf fn23 (* fn13 fn13))
  639.       (setf rfn13 (/ 1.0 fn13))
  640.       (setf w2r (+ (- coner (* zbr zbr)) (* zbi zbi)))
  641.       (setf w2i (- conei (* zbr zbi) (* zbr zbi)))
  642.       (setf aw2 (zabs w2r w2i))
  643.       (if (> aw2 0.25) (go label130))
  644.       (setf k 1)
  645.       (f2cl-lib:fset (f2cl-lib:fref pr (1) ((1 30))) coner)
  646.       (f2cl-lib:fset (f2cl-lib:fref pi_ (1) ((1 30))) conei)
  647.       (setf sumar (f2cl-lib:fref gama (1) ((1 30))))
  648.       (setf sumai zeroi)
  649.       (f2cl-lib:fset (f2cl-lib:fref ap (1) ((1 30))) 1.0)
  650.       (if (< aw2 tol) (go label20))
  651.       (f2cl-lib:fdo (k 2 (f2cl-lib:int-add k 1))
  652.                     ((> k 30) nil)
  653.         (tagbody
  654.           (f2cl-lib:fset (f2cl-lib:fref pr (k) ((1 30)))
  655.                          (-
  656.                           (*
  657.                            (f2cl-lib:fref pr ((f2cl-lib:int-sub k 1)) ((1 30)))
  658.                            w2r)
  659.                           (*
  660.                            (f2cl-lib:fref pi_
  661.                                           ((f2cl-lib:int-sub k 1))
  662.                                           ((1 30)))
  663.                            w2i)))
  664.           (f2cl-lib:fset (f2cl-lib:fref pi_ (k) ((1 30)))
  665.                          (+
  666.                           (*
  667.                            (f2cl-lib:fref pr ((f2cl-lib:int-sub k 1)) ((1 30)))
  668.                            w2i)
  669.                           (*
  670.                            (f2cl-lib:fref pi_
  671.                                           ((f2cl-lib:int-sub k 1))
  672.                                           ((1 30)))
  673.                            w2r)))
  674.           (setf sumar
  675.                   (+ sumar
  676.                      (* (f2cl-lib:fref pr (k) ((1 30)))
  677.                         (f2cl-lib:fref gama (k) ((1 30))))))
  678.           (setf sumai
  679.                   (+ sumai
  680.                      (* (f2cl-lib:fref pi_ (k) ((1 30)))
  681.                         (f2cl-lib:fref gama (k) ((1 30))))))
  682.           (f2cl-lib:fset (f2cl-lib:fref ap (k) ((1 30)))
  683.                          (*
  684.                           (f2cl-lib:fref ap ((f2cl-lib:int-sub k 1)) ((1 30)))
  685.                           aw2))
  686.           (if (< (f2cl-lib:fref ap (k) ((1 30))) tol) (go label20))
  687.          label10))
  688.       (setf k 30)
  689.      label20
  690.       (setf kmax k)
  691.       (setf zetar (- (* w2r sumar) (* w2i sumai)))
  692.       (setf zetai (+ (* w2r sumai) (* w2i sumar)))
  693.       (setf argr (* zetar fn23))
  694.       (setf argi (* zetai fn23))
  695.       (multiple-value-bind
  696.           (var-0 var-1 var-2 var-3)
  697.           (zsqrt sumar sumai zar zai)
  698.         (declare (ignore var-0 var-1))
  699.         (setf zar var-2)
  700.         (setf zai var-3))
  701.       (multiple-value-bind
  702.           (var-0 var-1 var-2 var-3)
  703.           (zsqrt w2r w2i str sti)
  704.         (declare (ignore var-0 var-1))
  705.         (setf str var-2)
  706.         (setf sti var-3))
  707.       (setf zeta2r (* str fnu))
  708.       (setf zeta2i (* sti fnu))
  709.       (setf str (+ coner (* ex2 (- (* zetar zar) (* zetai zai)))))
  710.       (setf sti (+ conei (* ex2 (+ (* zetar zai) (* zetai zar)))))
  711.       (setf zeta1r (- (* str zeta2r) (* sti zeta2i)))
  712.       (setf zeta1i (+ (* str zeta2i) (* sti zeta2r)))
  713.       (setf zar (+ zar zar))
  714.       (setf zai (+ zai zai))
  715.       (multiple-value-bind
  716.           (var-0 var-1 var-2 var-3)
  717.           (zsqrt zar zai str sti)
  718.         (declare (ignore var-0 var-1))
  719.         (setf str var-2)
  720.         (setf sti var-3))
  721.       (setf phir (* str rfn13))
  722.       (setf phii (* sti rfn13))
  723.       (if (= ipmtr 1) (go label120))
  724.       (setf sumbr zeror)
  725.       (setf sumbi zeroi)
  726.       (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  727.                     ((> k kmax) nil)
  728.         (tagbody
  729.           (setf sumbr
  730.                   (+ sumbr
  731.                      (* (f2cl-lib:fref pr (k) ((1 30)))
  732.                         (f2cl-lib:fref beta (k) ((1 210))))))
  733.           (setf sumbi
  734.                   (+ sumbi
  735.                      (* (f2cl-lib:fref pi_ (k) ((1 30)))
  736.                         (f2cl-lib:fref beta (k) ((1 210))))))
  737.          label30))
  738.       (setf asumr zeror)
  739.       (setf asumi zeroi)
  740.       (setf bsumr sumbr)
  741.       (setf bsumi sumbi)
  742.       (setf l1 0)
  743.       (setf l2 30)
  744.       (setf btol (* tol (+ (abs bsumr) (abs bsumi))))
  745.       (setf atol tol)
  746.       (setf pp 1.0)
  747.       (setf ias 0)
  748.       (setf ibs 0)
  749.       (if (< rfnu2 tol) (go label110))
  750.       (f2cl-lib:fdo (is 2 (f2cl-lib:int-add is 1))
  751.                     ((> is 7) nil)
  752.         (tagbody
  753.           (setf atol (/ atol rfnu2))
  754.           (setf pp (* pp rfnu2))
  755.           (if (= ias 1) (go label60))
  756.           (setf sumar zeror)
  757.           (setf sumai zeroi)
  758.           (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  759.                         ((> k kmax) nil)
  760.             (tagbody
  761.               (setf m (f2cl-lib:int-add l1 k))
  762.               (setf sumar
  763.                       (+ sumar
  764.                          (* (f2cl-lib:fref pr (k) ((1 30)))
  765.                             (f2cl-lib:fref alfa (m) ((1 180))))))
  766.               (setf sumai
  767.                       (+ sumai
  768.                          (* (f2cl-lib:fref pi_ (k) ((1 30)))
  769.                             (f2cl-lib:fref alfa (m) ((1 180))))))
  770.               (if (< (f2cl-lib:fref ap (k) ((1 30))) atol) (go label50))
  771.              label40))
  772.          label50
  773.           (setf asumr (+ asumr (* sumar pp)))
  774.           (setf asumi (+ asumi (* sumai pp)))
  775.           (if (< pp tol) (setf ias 1))
  776.          label60
  777.           (if (= ibs 1) (go label90))
  778.           (setf sumbr zeror)
  779.           (setf sumbi zeroi)
  780.           (f2cl-lib:fdo (k 1 (f2cl-lib:int-add k 1))
  781.                         ((> k kmax) nil)
  782.             (tagbody
  783.               (setf m (f2cl-lib:int-add l2 k))
  784.               (setf sumbr
  785.                       (+ sumbr
  786.                          (* (f2cl-lib:fref pr (k) ((1 30)))
  787.                             (f2cl-lib:fref beta (m) ((1 210))))))
  788.               (setf sumbi
  789.                       (+ sumbi
  790.                          (* (f2cl-lib:fref pi_ (k) ((1 30)))
  791.                             (f2cl-lib:fref beta (m) ((1 210))))))
  792.               (if (< (f2cl-lib:fref ap (k) ((1 30))) atol) (go label80))
  793.              label70))
  794.          label80
  795.           (setf bsumr (+ bsumr (* sumbr pp)))
  796.           (setf bsumi (+ bsumi (* sumbi pp)))
  797.           (if (< pp btol) (setf ibs 1))
  798.          label90
  799.           (if (and (= ias 1) (= ibs 1)) (go label110))
  800.           (setf l1 (f2cl-lib:int-add l1 30))
  801.           (setf l2 (f2cl-lib:int-add l2 30))
  802.          label100))
  803.      label110
  804.       (setf asumr (+ asumr coner))
  805.       (setf pp (* rfnu rfn13))
  806.       (setf bsumr (* bsumr pp))
  807.       (setf bsumi (* bsumi pp))
  808.      label120
  809.       (go end_label)
  810.      label130
  811.       (multiple-value-bind
  812.           (var-0 var-1 var-2 var-3)
  813.           (zsqrt w2r w2i wr wi)
  814.         (declare (ignore var-0 var-1))
  815.         (setf wr var-2)
  816.         (setf wi var-3))
  817.       (if (< wr 0.0) (setf wr 0.0))
  818.       (if (< wi 0.0) (setf wi 0.0))
  819.       (setf str (+ coner wr))
  820.       (setf sti wi)
  821.       (multiple-value-bind
  822.           (var-0 var-1 var-2 var-3 var-4 var-5)
  823.           (zdiv str sti zbr zbi zar zai)
  824.         (declare (ignore var-0 var-1 var-2 var-3))
  825.         (setf zar var-4)
  826.         (setf zai var-5))
  827.       (multiple-value-bind
  828.           (var-0 var-1 var-2 var-3 var-4)
  829.           (zlog zar zai zcr zci idum)
  830.         (declare (ignore var-0 var-1))
  831.         (setf zcr var-2)
  832.         (setf zci var-3)
  833.         (setf idum var-4))
  834.       (if (< zci 0.0) (setf zci 0.0))
  835.       (if (> zci hpi) (setf zci hpi))
  836.       (if (< zcr 0.0) (setf zcr 0.0))
  837.       (setf zthr (* (- zcr wr) 1.5))
  838.       (setf zthi (* (- zci wi) 1.5))
  839.       (setf zeta1r (* zcr fnu))
  840.       (setf zeta1i (* zci fnu))
  841.       (setf zeta2r (* wr fnu))
  842.       (setf zeta2i (* wi fnu))
  843.       (setf azth (zabs zthr zthi))
  844.       (setf ang thpi)
  845.       (if (and (>= zthr 0.0) (< zthi 0.0)) (go label140))
  846.       (setf ang hpi)
  847.       (if (= zthr 0.0) (go label140))
  848.       (setf ang (f2cl-lib:datan (/ zthi zthr)))
  849.       (if (< zthr 0.0) (setf ang (+ ang gpi)))
  850.      label140
  851.       (setf pp (expt azth ex2))
  852.       (setf ang (* ang ex2))
  853.       (setf zetar (* pp (cos ang)))
  854.       (setf zetai (* pp (sin ang)))
  855.       (if (< zetai 0.0) (setf zetai 0.0))
  856.       (setf argr (* zetar fn23))
  857.       (setf argi (* zetai fn23))
  858.       (multiple-value-bind
  859.           (var-0 var-1 var-2 var-3 var-4 var-5)
  860.           (zdiv zthr zthi zetar zetai rtztr rtzti)
  861.         (declare (ignore var-0 var-1 var-2 var-3))
  862.         (setf rtztr var-4)
  863.         (setf rtzti var-5))
  864.       (multiple-value-bind
  865.           (var-0 var-1 var-2 var-3 var-4 var-5)
  866.           (zdiv rtztr rtzti wr wi zar zai)
  867.         (declare (ignore var-0 var-1 var-2 var-3))
  868.         (setf zar var-4)
  869.         (setf zai var-5))
  870.       (setf tzar (+ zar zar))
  871.       (setf tzai (+ zai zai))
  872.       (multiple-value-bind
  873.           (var-0 var-1 var-2 var-3)
  874.           (zsqrt tzar tzai str sti)
  875.         (declare (ignore var-0 var-1))
  876.         (setf str var-2)
  877.         (setf sti var-3))
  878.       (setf phir (* str rfn13))
  879.       (setf phii (* sti rfn13))
  880.       (if (= ipmtr 1) (go label120))
  881.       (setf raw (/ 1.0 (f2cl-lib:fsqrt aw2)))
  882.       (setf str (* wr raw))
  883.       (setf sti (* (- wi) raw))
  884.       (setf tfnr (* str rfnu raw))
  885.       (setf tfni (* sti rfnu raw))
  886.       (setf razth (/ 1.0 azth))
  887.       (setf str (* zthr razth))
  888.       (setf sti (* (- zthi) razth))
  889.       (setf rzthr (* str razth rfnu))
  890.       (setf rzthi (* sti razth rfnu))
  891.       (setf zcr (* rzthr (f2cl-lib:fref ar (2) ((1 14)))))
  892.       (setf zci (* rzthi (f2cl-lib:fref ar (2) ((1 14)))))
  893.       (setf raw2 (/ 1.0 aw2))
  894.       (setf str (* w2r raw2))
  895.       (setf sti (* (- w2i) raw2))
  896.       (setf t2r (* str raw2))
  897.       (setf t2i (* sti raw2))
  898.       (setf str
  899.               (+ (* t2r (f2cl-lib:fref c (2) ((1 105))))
  900.                  (f2cl-lib:fref c (3) ((1 105)))))
  901.       (setf sti (* t2i (f2cl-lib:fref c (2) ((1 105)))))
  902.       (f2cl-lib:fset (f2cl-lib:fref upr (2) ((1 14)))
  903.                      (- (* str tfnr) (* sti tfni)))
  904.       (f2cl-lib:fset (f2cl-lib:fref upi (2) ((1 14)))
  905.                      (+ (* str tfni) (* sti tfnr)))
  906.       (setf bsumr (+ (f2cl-lib:fref upr (2) ((1 14))) zcr))
  907.       (setf bsumi (+ (f2cl-lib:fref upi (2) ((1 14))) zci))
  908.       (setf asumr zeror)
  909.       (setf asumi zeroi)
  910.       (if (< rfnu tol) (go label220))
  911.       (setf przthr rzthr)
  912.       (setf przthi rzthi)
  913.       (setf ptfnr tfnr)
  914.       (setf ptfni tfni)
  915.       (f2cl-lib:fset (f2cl-lib:fref upr (1) ((1 14))) coner)
  916.       (f2cl-lib:fset (f2cl-lib:fref upi (1) ((1 14))) conei)
  917.       (setf pp 1.0)
  918.       (setf btol (* tol (+ (abs bsumr) (abs bsumi))))
  919.       (setf ks 0)
  920.       (setf kp1 2)
  921.       (setf l 3)
  922.       (setf ias 0)
  923.       (setf ibs 0)
  924.       (f2cl-lib:fdo (lr 2 (f2cl-lib:int-add lr 2))
  925.                     ((> lr 12) nil)
  926.         (tagbody
  927.           (setf lrp1 (f2cl-lib:int-add lr 1))
  928.           (f2cl-lib:fdo (k lr (f2cl-lib:int-add k 1))
  929.                         ((> k lrp1) nil)
  930.             (tagbody
  931.               (setf ks (f2cl-lib:int-add ks 1))
  932.               (setf kp1 (f2cl-lib:int-add kp1 1))
  933.               (setf l (f2cl-lib:int-add l 1))
  934.               (setf zar (f2cl-lib:fref c (l) ((1 105))))
  935.               (setf zai zeroi)
  936.               (f2cl-lib:fdo (j 2 (f2cl-lib:int-add j 1))
  937.                             ((> j kp1) nil)
  938.                 (tagbody
  939.                   (setf l (f2cl-lib:int-add l 1))
  940.                   (setf str
  941.                           (+ (- (* zar t2r) (* t2i zai))
  942.                              (f2cl-lib:fref c (l) ((1 105)))))
  943.                   (setf zai (+ (* zar t2i) (* zai t2r)))
  944.                   (setf zar str)
  945.                  label150))
  946.               (setf str (- (* ptfnr tfnr) (* ptfni tfni)))
  947.               (setf ptfni (+ (* ptfnr tfni) (* ptfni tfnr)))
  948.               (setf ptfnr str)
  949.               (f2cl-lib:fset (f2cl-lib:fref upr (kp1) ((1 14)))
  950.                              (- (* ptfnr zar) (* ptfni zai)))
  951.               (f2cl-lib:fset (f2cl-lib:fref upi (kp1) ((1 14)))
  952.                              (+ (* ptfni zar) (* ptfnr zai)))
  953.               (f2cl-lib:fset (f2cl-lib:fref crr (ks) ((1 14)))
  954.                              (* przthr
  955.                                 (f2cl-lib:fref br
  956.                                                ((f2cl-lib:int-add ks 1))
  957.                                                ((1 14)))))
  958.               (f2cl-lib:fset (f2cl-lib:fref cri (ks) ((1 14)))
  959.                              (* przthi
  960.                                 (f2cl-lib:fref br
  961.                                                ((f2cl-lib:int-add ks 1))
  962.                                                ((1 14)))))
  963.               (setf str (- (* przthr rzthr) (* przthi rzthi)))
  964.               (setf przthi (+ (* przthr rzthi) (* przthi rzthr)))
  965.               (setf przthr str)
  966.               (f2cl-lib:fset (f2cl-lib:fref drr (ks) ((1 14)))
  967.                              (* przthr
  968.                                 (f2cl-lib:fref ar
  969.                                                ((f2cl-lib:int-add ks 2))
  970.                                                ((1 14)))))
  971.               (f2cl-lib:fset (f2cl-lib:fref dri (ks) ((1 14)))
  972.                              (* przthi
  973.                                 (f2cl-lib:fref ar
  974.                                                ((f2cl-lib:int-add ks 2))
  975.                                                ((1 14)))))
  976.              label160))
  977.           (setf pp (* pp rfnu2))
  978.           (if (= ias 1) (go label180))
  979.           (setf sumar (f2cl-lib:fref upr (lrp1) ((1 14))))
  980.           (setf sumai (f2cl-lib:fref upi (lrp1) ((1 14))))
  981.           (setf ju lrp1)
  982.           (f2cl-lib:fdo (jr 1 (f2cl-lib:int-add jr 1))
  983.                         ((> jr lr) nil)
  984.             (tagbody
  985.               (setf ju (f2cl-lib:int-sub ju 1))
  986.               (setf sumar
  987.                       (-
  988.                        (+ sumar
  989.                           (* (f2cl-lib:fref crr (jr) ((1 14)))
  990.                              (f2cl-lib:fref upr (ju) ((1 14)))))
  991.                        (* (f2cl-lib:fref cri (jr) ((1 14)))
  992.                           (f2cl-lib:fref upi (ju) ((1 14))))))
  993.               (setf sumai
  994.                       (+ sumai
  995.                          (* (f2cl-lib:fref crr (jr) ((1 14)))
  996.                             (f2cl-lib:fref upi (ju) ((1 14))))
  997.                          (* (f2cl-lib:fref cri (jr) ((1 14)))
  998.                             (f2cl-lib:fref upr (ju) ((1 14))))))
  999.              label170))
  1000.           (setf asumr (+ asumr sumar))
  1001.           (setf asumi (+ asumi sumai))
  1002.           (setf test (coerce (+ (abs sumar) (abs sumai)) 'double-float))
  1003.           (if (and (< pp tol) (< test tol)) (setf ias 1))
  1004.          label180
  1005.           (if (= ibs 1) (go label200))
  1006.           (setf sumbr
  1007.                   (-
  1008.                    (+ (f2cl-lib:fref upr ((f2cl-lib:int-add lr 2)) ((1 14)))
  1009.                       (* (f2cl-lib:fref upr (lrp1) ((1 14))) zcr))
  1010.                    (* (f2cl-lib:fref upi (lrp1) ((1 14))) zci)))
  1011.           (setf sumbi
  1012.                   (+ (f2cl-lib:fref upi ((f2cl-lib:int-add lr 2)) ((1 14)))
  1013.                      (* (f2cl-lib:fref upr (lrp1) ((1 14))) zci)
  1014.                      (* (f2cl-lib:fref upi (lrp1) ((1 14))) zcr)))
  1015.           (setf ju lrp1)
  1016.           (f2cl-lib:fdo (jr 1 (f2cl-lib:int-add jr 1))
  1017.                         ((> jr lr) nil)
  1018.             (tagbody
  1019.               (setf ju (f2cl-lib:int-sub ju 1))
  1020.               (setf sumbr
  1021.                       (-
  1022.                        (+ sumbr
  1023.                           (* (f2cl-lib:fref drr (jr) ((1 14)))
  1024.                              (f2cl-lib:fref upr (ju) ((1 14)))))
  1025.                        (* (f2cl-lib:fref dri (jr) ((1 14)))
  1026.                           (f2cl-lib:fref upi (ju) ((1 14))))))
  1027.               (setf sumbi
  1028.                       (+ sumbi
  1029.                          (* (f2cl-lib:fref drr (jr) ((1 14)))
  1030.                             (f2cl-lib:fref upi (ju) ((1 14))))
  1031.                          (* (f2cl-lib:fref dri (jr) ((1 14)))
  1032.                             (f2cl-lib:fref upr (ju) ((1 14))))))
  1033.              label190))
  1034.           (setf bsumr (+ bsumr sumbr))
  1035.           (setf bsumi (+ bsumi sumbi))
  1036.           (setf test (coerce (+ (abs sumbr) (abs sumbi)) 'double-float))
  1037.           (if (and (< pp btol) (< test btol)) (setf ibs 1))
  1038.          label200
  1039.           (if (and (= ias 1) (= ibs 1)) (go label220))
  1040.          label210))
  1041.      label220
  1042.       (setf asumr (+ asumr coner))
  1043.       (setf str (* (- bsumr) rfn13))
  1044.       (setf sti (* (- bsumi) rfn13))
  1045.       (multiple-value-bind
  1046.           (var-0 var-1 var-2 var-3 var-4 var-5)
  1047.           (zdiv str sti rtztr rtzti bsumr bsumi)
  1048.         (declare (ignore var-0 var-1 var-2 var-3))
  1049.         (setf bsumr var-4)
  1050.         (setf bsumi var-5))
  1051.       (go label120)
  1052.      end_label
  1053.       (return
  1054.        (values nil
  1055.                nil
  1056.                nil
  1057.                nil
  1058.                nil
  1059.                phir
  1060.                phii
  1061.                argr
  1062.                argi
  1063.                zeta1r
  1064.                zeta1i
  1065.                zeta2r
  1066.                zeta2i
  1067.                asumr
  1068.                asumi
  1069.                bsumr
  1070.                bsumi)))))
  1071.  
  1072.